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Abstract 

Models of continuous opinion dynamics under bounded confidence 
have been presented independently by Krause and Hegselmann and by 
Deffuant et al in 2000. They have raised a fair amount of attention in 
the communities of social simulation, sociophysics and complexity sci- 
ence. The researchers working on it come from disciplines as physics, 
mathematics, computer science, social psychology and philosophy. 

In these models agents hold continuous opinions which they can grad- 
ually adjust if they hear the opinions of others. The idea of bounded 
confidence is that agents only interact if they are close in opinion to each 
other. Usually, the models are analyzed with agent-based simulations in 
a Monte-Carlo style, but they can also be reformulated on the agent's 
density in the opinion space in a master-equation style. The contribu- 
tion of this survey is fourfold. First, it will present the agent-based and 
density-based modeling frameworks including the cases of multidimen- 
sional opinions and heterogeneous bounds of confidence. Second, it will 
give the bifurcation diagrams of cluster configuration in the homogeneous 
model with uniformly distributed initial opinions. Third, it will review the 
several extensions and the evolving phenomena which have been studied 
so far, and fourth it will state some open questions. 

opinion dynamics, continuous opinions, cluster formation, bifurcation patterns 
PACS Nos.: 89.20.-a; 89.65.-s 

1 Introduction 

The term "opinion dynamics" nowadays summarizes a wide class of different 
models differing in heuristics, formalization as well as in the phenomena of 
interest. The latter range from emergence of fads, minority opinion spreading, 
collective decision making, finding and not finding of consensus, emergence of 
political parties, minority opinion survival, emergence of extremism and so on. 
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This paper deals with these phenomena in models of continuous opinion 
dynamics under bounded confidence. 'Continuous' refers to the opinion issue 
and not to the time. Thus opinions in continuous opinion dynamics should be 
expressable in real numbers where compromising in the middle is always possi- 
ble. Example issues are prices, tax rates or predictions about macroeconomic 
variables. The political spectrum is also often mapped to a continuum from left 
to right wing. 

In opinion dynamics one considers a set of agents where each holds an opinion 
from a certain opinion space. She may change her opinion when she gets aware 
of the opinions of others. In the physics literature discrete opinion spaces (clas- 
sically binary opinions) (JJ 0; d) have dominated research due to their striking 
analogy with spin systems. Sometimes they have been extended to more than 
two spin values, which are ordered and thus get closer to continuous opinion 
dynamics fill). 

In recent years two models of genuinely continuous opinion dynamics under 
bounded confidence have raised the interest of the sociophysics community: the 
models of Hegselmann and Krause (0; 0; H) and Deffuant, Weisbuch and others 
(9; 10) . Due to physicist's research some progress in understanding the dynamics 
of these models has been made, especially by introducing the master equation 
on the agents density in the opinion space for these type of models. For a recent 
review in this broader context of opinion dynamics including discrete opinions 
from a physicist's perspective see (11). 

Let us consider a population of agents which hold diverse opinions about 
certain issues expressible in real numbers. Each agent is willing to change her 
opinion if she hears the opinions of others by adjusting towards those opinions. 
Every adjustment in terms of averaging is possible due to the continuous nature 
of the opinions. Further on, consider agents to have bounded confidence. That 
means an agent is only willing to take those opinions into account, which differ 
less than a certain bound of confidence e from her own opinion. 

The Deffuant- Weisbuch (DW) model and the Hegselmann-Krause (HK) model 
both rely on the idea of repeated averaging under bounded confidence. They 
differ in their communication regime. In the DW model agents meet in random 
pairwise encounters after which they compromise or not. In the HK model, each 
agent moves to the average opinion of all agents which lie in her area of confi- 
dence (including herself). Actually, the DW model contains another parameter 
which controlls how close an agent moves to the opinion of the other. But it 
has turned out that this parameter has only an effect on convergence time in 
the basic model. Therefore we neglect it in the basic analysis and discuss its 
impact afterwards. 

The DW model was partly inspired by the famous Axelrod model about the 
dissemination of culture (|12l ) where something similar to the bounded confidence 
assumption is implemented. It has been developed in a project about improving 
agri-environmental policies in the European union. 

The HK model has been presented by Krause 0) in a mathematical 
context as a nonlinear version of older consensus models \l4k ll5T ). It was 
analyzed through computer simulations by Hegselmann and Krause (|8() in the 
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context of social simulation and has gained a lot of attention since. 

The next section presents the two models in their original agent-based ver- 
sions and in their density-based formulation inspired by statistical physics. Sec- 
tion 3 shows and explains the bifurcation diagrams of both models which serve 
as a reference for the review of several extensions in Section 4. 

2 The models 

From the heuristic description of the models one can either define agent-based 
dynamics for a finite population of n agents or density-based dynamics for a 
density function which determines the agents density in the opinion space. The 
latter approach can then be interpreted as taking the infinite limit n — > oo of 
agent-based models. It is a classical tool of statistical physics often called the 
derivation of a master equation or rate equation which was first applied to social 
systems bei Weidlich (fl(lli"7h. 

Nevertheless, the basic ingredient of both approaches is a continuous opinion 
space S C Mr. The d real numbers represent opinions on d different subjects. 
Usually, only compact and convex sets are regarded as appropriate opinion 
spaces. The simplest case that we will discuss is the one-dimensional interval 
[0,1]. But definitions extend naturally to more dimensions. 

In an agent-based model the state variable is an opinion profile x(t) G S n , 
which is a vector of vectors. For agent i, it contains the opinion vector x % G S. 
For an initial opinion profile x(0) dynamics are defined recursively as x(t + 1) = 
f(t,x(t)). The pair (S, f) is then a discrete dynamical system. 

In a density based model, the state variable is a density function on the 
opinion space P(t,-) : S — » R> with J s P(t,x)dx = 1 for all t. Dynam- 
ics are then defined for a given initial density function P(0, •) as the evolu- 
tion of the density function in time. Time is sometimes regarded as discrete 
and sometimes as continuous. In the continuous-time case we consider the 
differential equation JjP(i, x) = g(P(t,-)) where g operates on the space of 
density functions. Then we try to find solutions of the differential equation 
analytically or with numerical solvers. In the discrete-time case we replace 
the differential operator by the difference operator and write AP(t + l,x) = 
P(t + l,x) — P(t,x) = g(P(t, •)). We can then directly compute the solution 
recursively P(t+ l.x) = P(t,x) + AP(t+l). Actually, both discrete {II|;[ig|) and 
continuous (|20l ; |2 lh time were used to compute solutions numerically. It turned 
out that discrete and continuous solutions of the DW model are quite similar 
but discrete and continuous solutions of the HK model differ qualitatively, which 
we will point out later. 

In the following we will define the agent-based DW and HK model in detail, 
as well as their density-based counterparts. The definitions will already include 
the extensions of multidimensional opinions and heterogeneous bounds of con- 
fidence. In a first step it is enough to think of the opinion space as the interval 
[0, 1] and to regard all bounds of confidence as equal. 
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2.1 The Deffuant-Weisbuch model 



We proceed with the definition of the agent-based version of the DW model. 

Definition 1 (Agent-based DW model) Let there be n £ N agents and an 
appropriate opinion space S C Mr. Given an initial profile x(0) £ S n , bounds 
of confidence ei,...,e„ > 0, and a norm \\-\\ we define the agent-based DW 
process as the random process (x(t))t^ that chooses in each time step t £ N 
two random agents i,j which perform the action 

^ = f (x»'(t)+i*(t))/2 <Si 

1 x l {t) otherwise. 



The same for x° (t + 1) with i and j interchanged. 

If Ex = ••• = £« we call the model homogeneous, otherwise heterogeneous. 

It has been shown (|22n that the homogeneous process always converges to 
a limit opinion profile (usually not in finite time). The same is observed in 
simulations of the heterogeneous case, but a proof is lacking. A limit profile in 
the homogeneous case is an opinion profile where for each two opinions Xi , Xj it 
holds that they are either equal (belong to the same cluster) or have a distance 
larger than e. Thus further changes are not possible regardless of the choice of 
i and j. Further on, it is easy to see that the average opinion over all agents is 
conserved during dynamics (23), but only in the homogeneous case. 

For example trajectories of the process see $t liH). 

For the density-based DW model we extend the model defined in (j20l ) to 
populations of agents with heterogeneous bounds of confidence. 



Definition 2 (Density-based DW model) Let S C 



be an appropriate 



is of 
1.0 



opinion space, [e\ , e 2 ] be an interval of possible bounds of confidence and the 
initial density function on the opinion space times the interval of bounds 
confidence be P(0, •, •) : S x [£1,22] — > [0,oo]. with J s J^ 2 dxdeP(Q, x, e) 
For abbreviation we define the aggregated density as P(t 7 x) — deP(t,x,e). 
With this we define the differential equation 

—P(t,x,e)= J dxi J de ( / dx 2 P(t,xi,e)P(t,x 2 ) + 



\\xi-x 2 \\<e 



P(t,xi)P(t,x 2 ,e) S(x 



X\ - x 2 



dx 2 (P(t, x 1 ,e)P(t, x 2 )S(x - x^ + P(t, xi)P(t, x 2 ,e)8(x - x 2 )) 



II a=l — a=2 II <^ 



(1) 



1 P(0, ■,■) is a density function at time zero over the opinion space and the interval of 
bounds of confidence, P(0, x, e)dxde represents the proportion of agents which hold opinions 
in [x, x + dx] and bounds of confidence in [e, e + de]. 
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The continuous-time density-based DW process is then defined as the solution 
of the initial value problem. The discrete-time density-based DW process is the 
sequence (P(t, ■, -))teN recursively defined by 

P{t + 1, x, e) = P(t, x, e) + AP(t, x, e) 

with AP(t,x,e) = j^P(t,x,e). 

Notice that it is important to distinguish the aggregated density Pit, x) 
defined above and the density for a given bound of confidence Pit, x, e) to un- 
derstand how an agent looks at all agents regardless of their bound of confidence 
and adjusts to the ones within her bound of confidence. 

For a homogeneous bound of confidence e — [e%, e 2 ] the dynamical Equation 
(fTJ) reduces to 

^P(t,x) = J dx x J dx 2 ( P(t, Xl )P(t,x 2 )( 25{x~ Xl+ 2 X2 ) 

s \\xi-w 2 \\< e >- — rr. — — ^ 

traction joining state x 

- {5{x - x{) + 5{x - x 2 )\ ) ) 

fraction leaving state x 

(2) 

which has been reported first in (20I ) and gives a good starting point for under- 
standing populations with heterogeneous bounds of confidence. The ^-functions 
represent the position where the mass of the two opinions x\ and x 2 jumps to. 
Equation ([2|) can be transformed to a form free of ^-distributions, 

^P(t,x)= J dy4P(t,y)P(t,2x-y) J dy2P(t,y). (3) 

\\x-y\\<§ \\x-y\\<s 

One insight is that the gain term for opinion a; is a 'bounded convolution' of 
Pit, •) with itself at the point 2x. 

It is pointed out in (|20r ) that Equation ^ conserves the mass and the mean 
opinion of the density. For Equation |T]) this holds only for the mass, which 
is conserved for each value of e and of course for the aggregated density, too. 
This gives rise to the fact that under heterogeneous bounds of confidence overall 
drifts (in terms of a drift of the mean opinion) to more extremal opinions may 
happen. This has been simulated by lj24f ) in a very stylized model. 

In ( 20) it is shown for S = [0, 1] and e > 0.5 that P(t, •) in Equation © 
goes with t — > 00 to a limit distribution P(oo, x) = <5(.t — xq) where Xq is the 
mean opinion. For lower epsilon it is postulated (and verified by simulation) 
that P(oo,x) = y^j—! mjSjx — Xi) where r is the number of evolving opinion 
clusters, Xi £ [0, 1] is the position of cluster i and rrii is its mass. All clusters 
must fulfill the conservation laws X)i=i m i — 1 an d x % m i is equal to the 

conserved mean opinion. Further on, all different clusters i ^= j must fulfill 
\xi — Xj\ > e. This conicides with the result for the agent-based model. 
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In general the right hand side of Equation @ does not depend continuously 
on P(t, •). El But this happens only in the case where P(t, •) contains (5-functions 
and is thus not a 'normal' function in x but a distribution. Noncontinuity of 
the right hand side of ||2J| implies that there is no or at least no unique solution. 
Fortunately, it turned out that ^-functions do not evolve during the process but 
only in the limit t — > oc. If we start with P(t, •) without ^-functions trajectories 
can be computed as it is done in |20h with a fourth order Adams-Bashforth 
algorithm. 

Moreover, it turned out in comparing results of (fiol ) and (flit ) that the 
discrete-time density-based process leads to nearly the same limit densities as 
the continuous-time density-based process. In (flit flit ), additional to discrete 
time, the space P(t, •) has been discretized to p(t) € R" . Dynamics can then 
be defined as an interactive Markov chain (a state and time discrete Markov 
chain where transition probabilities depend on the actual state). Essentially, 
the interactive Markov chain is the same what evolves by discretization of the 
continuous opinion space for numerical computation in (|20 ). 

In a nutshell, there are different approaches but results for the density-based 
DW model lie all close together. This does not hold for the HK model as we 
will see. 

Examples for trajectories of the process can be found in ( TiJ 20l ). 



2.2 The Hegselmann-Krause model 

Here we will give the HK model, first in its original agent-based version. 

Definition 3 Agent-based HK model Let there fee n S N agents and an 
appropriate opinion space S C K d . 

Given an initial profile x(0) € S n , bounds of confidence e\, . . . ,e n > and 
a norm \\-\\ we define the HK process (x{t)) t ^n recursively through 

x(t + 1) = A{x(t), ei, ... , s n )x(t), (4) 

with A(x,Si, . . . , e„) being the confidence matrix defined 

a ( \ / *i\ x) $3 G Iei(i,x) 

A ij (x,e 1 ,...,e n ):=i * 1 ^ x > 

I U otherwise, 

with I £i (i,x) := {j S n | \\x l — x 3 : || < e^} being the confidence set of agent i with 
respect to opinion profile x. In other words, the new opinion x l {t + 1) is the 
arithmetic average over all opinions in x(t) differing from it by not more than 

If £\ = ■ ■ ■ — e n we call the model homogeneous, otherwise heterogeneous. 



2 Consider P(t, ■) = — 1) + ^S(x + 1 + f) then we have a fixed point for e = 2 

and any £ > 0. But for § = applying the right hand side of Equation JJJ would give 
\8{x - 1) + |<5(x) + §5(3 + 1). 
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It has been shown (|25l ; |26( ) that the homogeneous process always converges 
to a limit opinion profile in finite time. The same convergence is observed in 
simulations for the heterogeneous case but it does not occur in finite time and a 
proof is lacking. A limit profile in the homogeneous case is always a fixed point 
x* = A(x*,e)x*. In x* it holds for each pair of two opinions x % ,x^ that they 
are either equal (belong to the same cluster) or have a distance larger than s. 
The proof for this is not difficult but also not trivial and can be found in (flil ). 
In contrast to the DW model the mean opinion is not conserved. The mean 
opinion is only conserved when the initial profile is symmetric around its mean 
opinion. 

For example trajectories of the process see (fil; liH). 

We now turn our attention to the density-based version of the HK model for 
which we extend models independently developed in |2lh and (flih . 

For the definition of the density-based HK model we need the definition of 
the e-local mean 

C*yP(t,y)dy 



M 1 (x,P(t,-),s) 



i:::p(t,y)dy 



It gives the expected value of P(t, •) in the interval [x — e, x + e]. The nominator 
gives the first moment in the e-interval around x while the denominator gives 
the necessary renormalization by the probability mass in that interval. 

Definition 4 (Density-based HK model) Let S C M d be an appropriate 
opinion space, [61,62] be an interval of possible bounds of confidence and the 
initial density function on the opinion space times the interval of bounds of 
confidence be P(0, •, •) : S x [£1,62] — > [0, 00] with J s dxdeP(0,x,e) = 1. □ 
For abbreviation we define the aggregated density as P(t,x) — deP(t 7 x,e). 
With this we define the differential equation 



d f 

—P(t,x,e) = dy 



£2 

deP(t,y,e)5(x - Mi(y,P(t,-),e)) 



£1 



P(t,x,e)S(y-M 1 (x,P(t,-),s)) (5) 



The continuous-time density-based HK process is then defined as the solution 
of the initial value problem. The discrete-time density-based HK process is the 
sequence (P(t, •, -))ten recursively defined by 

P{t + 1, x, e) = P(t, x, e) + AP(t, x, e) 

with AP(t,x,e) = ^P(t,x,E). 



3 P(0, ■,■) is a density function at time zero over the opinion space and the interval of 
bounds of confidence, P(0, x, e)dxde represents the proportion of agents which hold opinions 
in [x,x + dx] and bounds of confidence in [e, e + de]. 
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As for the density-based DW model it is important to distinguish the ag- 
gregated density P(t,x) defined above and the density for a given bound of 
confidence P(t, x, e) to understand how an agent looks at all agents regardless 
of their bound of confidence and takes into account all within her bound of 
confidence. 

For a homogeneous bound of confidence e = [ei, £2] the dynamical Equation 
||SJ) reduces to 



° P(t, x) = I 5(Mi(y, P(t, •),£)- x)P(t, y) - SiM^x, P(t, -),e)- y)P(t, x) dy. 

fraction joining state x fraction leaving state x 

(6) 



dt 



This can be transformed to the form 



J dyS(y)- 



P(t,y) 



\M>(y,P(t,-),e)\ 

\{y is root of x-M 1 (y,P(t,-),s)} 



P(t,x) 



(7) 



with the M[ being the derivative with respect to y. The function M\{y, P(t, •), e) 
is monotone in y for all P. Therefore x — Mi(y, P(t, -),s) has usually only one 
or no root. In the case of one root the integral in Equation ([7]) reduces to 
with yo being that root. But there might be a continuum of 



P(yo) 



\M[(y ,P(t,-),e)\ 

roots due to a plateau of x — Mi(y,P(t, ■),£) at zero, because M\(y,P(t, -),e) 
is not strongly monotone. 

It is shown in (|2ll ) that Equation © conserves the mass. Further on it is 
claimed that the mean opinion is also conserved. Although the derivation in 
mi) relies on the symmetry of the initial opinion density around its mean the 
authors forgot to mention that the mean opinion is not generally conserved for 
other initial conditions. 

In general the right hand side of Equation (J5]) does not depend continuously 
on P(t, •). (See example above.) But unfortunately a 5-peak in the density 
may evolve in this model by computing the right hand side of Equation ([6]) 
even with a given P(t, ■) which does not contain a (5-peak. Therefore the 
existence of unique solutions to j6| is not generally assured for every initial 
condition. Nevertheless they were computed in |2lh with a fourth order Runge- 
Kutta algorithm and smooth looking dynamics were reached. 

But the picture changes drastically when we switch to the discrete-time 



density-based HK model, which has been studied in (|18l : l27t Il9f ) as an interactive 
Markov chain with 1000 and 1001 opinion classes. This coincides with the 
discretization of the interval [0, 1] into 1000 bins in (I2l[). In the discrete-time 
case two phenomena occur which were not reported in (|2 lh : Consensus striking 



4 Consider P(t, •) to be uniformly distributed on [0, U 1] and e = 0.5. Then all the 
'agents' with opinions in the interval [0, have e-local mean equal to g at that position we 
will have a <5-peak. 
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back for lower e after a phase of polarization (1271) and the importance of an 
odd versus an even number of opinion classes (|18(). Under an odd number, 
consensus is much easier to reach because under an even number, the central 
mass contracts into two bins which might be divided more easily compared to 
the odd case. In the odd case the mass is stored in one bin, and mass in one bin 
can only move jointly by definition of the process. We discuss the differences of 
continuous-time and discrete-time results in the next section. 

A numerical issue which is not clearly defined in (|21l ) is how moving mass is 
assigned to bins when computing the right hand side of Equation ©. Usually, 
the e-local mean of a given bin lies not directly in a bin but between two bins. 
So, there has to be a rule how to distribute it between these two bins. In 
(flU 27; [l9|) it is done proportionally to the distances from each bin. If the 



e-local mean lies closer to one bin, this bin receives more of the moving mass. 

In this difference of continuous versus discrete time, the HK model distin- 
guishes from the DW model. Nevertheless, the possible fixed points of the 
dynamics are the same as in the DW model. It is again every distribution like 
P(oo, x) = Yli=i Triid(x — Xi) with \xi — Xj\ > e for all i ^ j and J2i=i m i = !• 
This coincides with the result for the agent-based model. But for both models 
it is not easy to determine the limit density out of the initial density. Nearly 
every study until now has studied only uniformly distributed initial densities. 



Example trajectories of the process can be found in (21; 1 



3 Bifurcation Diagrams 

In this section we give the basic bifurcation diagrams for both homogeneous 
density-based models. A bifurcation diagram shows the location of clusters in 
the limit density versus the continuum of values of the bound of confidence 
e. So, one can determine the attractive cluster patterns for each bound of 
confidence and observe transitions of attractive patterns at critical values of e. 
The bifurcation diagrams serve as references for discussions on robustness to 
extensions of the homogeneous models and the impact of these extensions, e.g. 
on the critical bounds of confidence. 

The top plot of Figure [T] shows the bifurcation diagram for the homogeneous 
density-based DW model with uniform initial density in the opinion space [0, 1]. 

A bifurcation diagram is made by computing the trajectory of the corre- 
sponding interactive Markov chain with 1001 opinion classes until clusters have 
evolved which spread not more than e and are further away from each other than 
s. Clusters are determined by collecting intervals in the opinion space where 
the density is positive. Therefore, one has to decide on a threshold for which 
mass is regarded as zero. This threshold was 10 -9 here. (For computational 
details see (fll ).) 

The bifurcation diagram shown here resembles the bifurcation diagram shown 
m (20) but with a transformation of variables as A = ^. In (j2Qj) A is the 
changing variable and determines the opinion space [—A, A] and the bound of 
confidence is fixed to one. Here we have a fixed opinion space and a changing 
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bound of confidence because it is always done that way in agent-based simula- 
tions and this framework can be extended to models with heterogeneous bounds 
of confidence. 

The fat lines in the top plot display the locations of clusters which have 
a major mass, the thin lines the location of clusters with a minor mass and 
the medium line the central cluster. Dotted lines are just for orientation, they 
show the interval [0.5 — e, 0.5 + e]. The two lower plots show the masses of the 
related clusters. The most striking feature is the existence of minor clusters 
at the extremes and between major clusters. In agent-based simulations these 
minor clusters often do not evolve because of low population sizes, nevertheless 
outliers are frequently reported. These outliers make it difficult to count the 
number of major clusters correctly. In [l(J authors therefore decided to not 
count clusters with only one agent. The bifurcation diagram in |20l ) clarified 
that minor clusters exist for structural reasons. 

For e > 0.5 only one big central cluster evolves. As e decreases bifurcations 
and nucleations of clusters occur. First, the nucleation of two minor extremal 
clusters, then the bifurcation of the central cluster into two major clusters, and 
third the rebirth of the central cluster. For details see (19|; 2(3). This bifurcation 
pattern is then repeated in shorter e-intervals. The length of these intervals 
seems to scale with K This can be better seen in the bifurcation diagram 
of ©, where the bifurcation pattern seems to repeat itself on intervals that 
converge towards a length of about 2.155. But this constant has only been 
derived numerically and it is not clear that the bifurcation pattern will repeat 
that regularly, although it looks very much like that. This result resembles the 
rough Tjj-rule reported in (0; [lfl ) for agent-based simulation, which says that 
the number of major clusters after cluster formation is roughly determined as 
the integer part of 

A certain e-phase of interest is e S [0.27, 0.35]. As e decreases the mass of the 
minor cluster decreases, too. Thus the mass of the central cluster must increase 
slightly until the central cluster bifurcates. Thus, surprisingly the central cluster 
has a larger mass very close to the critical e. For further minor clusters the 
same intermediate low-mass region exist. It even looks as if the minor cluster 
disappears for a short intermediate e-range. The question if these gaps exist 
has also been posed in (|20l ). 

The top plot of Figure [2] shows the bifurcation diagram for the homogeneous 
density-based HK model with uniform initial density in the opinion space [0, 1]. 
The bifurcation diagram was computed the same way as for the DW model. 
But the detection of the final density was much easier because the time-discrete 
process converges in finite time. (For computational details see (HH).) 

The lines in the top plot display the locations of clusters. Unlike the DW 
model the HK model has no minor clusters. The lower plot shows the masses of 
the related clusters. The convergence time (which is finite here) is the subject 
of the lowest plot. There are certain e- values were convergence time is really 
long. The plot shows a zoom, which does not display the longest times. The 
longest time was 21431. 

For about e > 0.19 consensus is reached, but for values only slightly above, 
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DW model: bifurcation diagram 
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Figure 1: Bifurcation diagram for the DW model 
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it takes a long time. (See (21) for a discussion why or (flil ) for explanatory 
figures of the processes.) So, e = 0.19 is the consensus threshold. In (f2lh the 
consensus threshold lies at about e — 0.22. This difference might have two 
reasons. First, it could be caused by the differences of the continuous-time and 
the discrete-time process, and second it could result from the fact that they 
choose 1000 bins, so an even number. The computation in Figure [H wa s done 
with 1001 bins. This may really make a difference as pointed out in (|18T). 

A striking feature of the bifurcation diagram is that after an e-phase of 
polarization into two clusters (and a tiny central cluster) there is a phase of even 
lower e where we reach consensus again at the expense of very long convergence 
time. Details about this phenomenon can be found in |27h. Again this feature 
does not evolve in the bifurcation diagram of (|2lf). This might have again two 
reasons. First, the continuous-time approach, and second their termination 
criteria was so rough that it neglects the small movements which bring the two 
intermediate clusters together. Indeed there is a small hint for this reason in 
the bifurcation diagram in |2ll): There is a small hill in the convergence time 
at the region where convergence to consensus happens in Figure [5] 

So, the HK bifurcation diagram differs a lot from the DW bifurcation di- 
agram. Nevertheless, they are similar in the fact that the bifurcation pattern 
seems to repeat on time scales scaling with ~. Again, there is no proof of this 
universal scaling. Moreover there is evidence that the regularity gets destroyed 
for lower e. (See (|27l ) for details, the evidence can be seen by studying onesided 
dynamics.) 

Figure [3] shows both bifurcation diagrams in one plot to compare the dif- 
ferences. (Black lines are DW model, red (color online) lines ar HK model.) 
Notice the different critical values of e for the consensus transition. The con- 
sensus transition here is the bifurcation point between two major clusters and a 
big central cluster. So, we neglect the minor clusters in the DW model for this 
definition. 

The DW and HK model have been studied with several extensions, which we 
will review in the following with respect to the bifurcation diagrams presented 
here. 



4 Extensions 

The Definitions 1, 2, 3 and 4 already contain the possibility of multidimensional 
opinions and heterogeneous bounds of confidence for agents, respectively sub- 
populations. Several further extensions which are based on continuous opinions 
and bounded confidence haven been proposed and analyzed to some extent. A 
good example of the overwhelming number of ideas which need to be checked 
is (j28l ). We will review the known results and propose some open questions. 

Nonuniform initial conditions This topic is not really an extension 
but a very important issue. Nearly every study used random and uniformly 
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HK and DW model: bifurcation diagram 
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Figure 3: Bifurcation diagram for the DW and HK model 



distributed initial opinion profiles or initial densities which are uniformly dis- 
tributed in the opinion space. The right classification of initial profiles leading 
to structurally different final cluster configurations remains a task for future 
work. At the agregate level the agent-based simulations are in fair agreement 
with the density-based processes, although random fluctuations play important 
roles in individual trajectories even for large populations. This topic is also 
discussed in (f29h . 



Multidimensional opinions Until now, we have only considered the 
opinion space [0, 1]. Other lengths of intervals will not deliver new insight, be- 
cause they will give the same structural behavior only on scaled e-levels. So, we 
have a complete picture in studying the unit interval. This changes when inves- 
tigating multidimensional opinions. The usual first step is to study S = [0, l] 2 , 
see (|2ll l26t 1301). Final clusters then form a regular square latice but with some 
interferences in the center, where a central cluster may or may not exist. See 
(30) for a circular disk as the opinion space in an agent-based HK model and 
how it can have a positive impact on finding consensus by preventing corner 
clusters. See (jicj ) for a hexagonal and a triangular opinion spaces. 

Another interesting parameter, which plays a role in more dimensions is 
the distances measure which is a norm of the difference of two opinions in this 
paper, but of course any metric might be a candidate. The most common p- 
norms (||a;|| = Ei=i I 2 '*!*')') have been studied. These are p = 1, 2, oo in 
(j26l ) for the agent-based HK model (see also |30l)) and p = 2, oo in (21) for the 
density-based HK model. There is some impact of the norm but the interplay of 
the distance measure and the shape of the opinion space is not well understood. 
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There are only two studies about the impact of the dimensionality parameter 



d itself. In |27() ci-dimensional unit simplexes (S = A d = {x e M>V | Ei=i 
1}) as opinion spaces were studied up to d = 7. The simplex opinion space repre- 
sents opinions which are under budget constraints. One opinion is a distribution 
of resources to d + 1 projects (the opinion issues) but with a fixed amount of 
resources. It has been shown that raising d lowers the consensus threshold in 
the agent-based DW as well as in the agent-based HK model. But the consensus 
threshold cannot be pushed down to zero. So, there is a limit above which rais- 
ing the number of issues d + 1 brings only marginal positive effects on fostering 
the chances for consensus. 

The budget constraints are crucial for this positive effect as it is shown in 



(311). If opinions are not under budget constraints like in S = [0, 1] raising d 
raises the consensus threshold. This has been checked for DW and HK commu- 
nication, 1- and oo-norm. 

In a nutshell, if you want to foster consensus, bring more interrelated issues 
into the discussion. If you want to prevent consensus, bring more independent 
issues into the discussion. 

One source of this effect is the growth of the volumes of opinion spaces com- 
pared with the growth of volumes of areas of confidence with rising dimension. 
Under budget constraints the volume of the area of confidence growth faster 
than the volume of the opinion space. It is the other way round without budget 
constraints. 

Another interesting new approach in multidimensional continuous opinion 
dynamics is to study dynamics of structured attitudes and opinions ([H), where 
attitudes are multidimensional but opinions are an aggregation of attitudes and 
dynamics happen on both levels. 

Heterogeneous bounds of confidence The first brief study of hetero- 
geneous bounds of confidence has been done in (fioh for the DW model with a 
population of agents with two levels of bounds of confidence, the open-minded 
and the closed minded agents. They studied a set of 192 closed-minded agents 
(si = 0.2) and eight open-minded agents (£2 = 0-4) and saw that in the short 
run, the cluster pattern of the closed-minded (two big clusters) dominated while 
in the long run, the one of the open-minded (consensus) evolved. 

Meanwhile it has been shown by analyzing the density-based DW model flil ) 
that there is even more impact of heterogeneity. Consensus can be achieved by 
mixing closed- and open-minded agents even if both bounds of confidence are 
far below the critical value of the consensus transition (e.g. £1 = 0.11 and 
£2 = 0.22). But on the other hand drifting of clusters of open-minded towards 
clusters of closed-minded agents becomes a generic feature of dynamics which 
can amplify some asymmetric disturbances in the initial profile. So in the end 
a final consensual cluster may lie very far from the initial average opinion. 

This drifting phenomenon is the source of the drift to the extremes studied in 
HiiS). There, very few closed-minded agents are initially set to the extremes 



of the opinion space and open-minded agents are set in between. So, very stable 
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clusters of extremists always evolve at the extremes and the question is where 
the central agents go. They may form a central cluster, divide to both extremes 
or drift together to one extreme, depending on the threshold of the open-minded 
agents. 

The HK model shows similar drifting phenomena when heterogeneous bounds 
of confidence are introduced. Further on, another phenomenon is 'sitting be- 
tween the chairs' of the open-minded agents. It might happen that two clusters 
of closed-minded agents form and open-minded agents form a cluster in between 
although they take the opinions of both closed- minded clusters into account. 



Social networks and communication regimes Another issue which 
has been studied is the influence of social networks. Communication is not 
only restricted by the bounded confidence assumption but also by a given static 
social network. In |34|) and (35) the agent-based DW and HK model has been 
set up on networks with scale-free link distribution a la Barabasi- Albert, random 
graphs a la Erdos-Renyi, and square lattices. They gave evidence for the fact 
that in the DW model for large enough networks consensus is always reached for 
e > 0.5, which was later proven in (|36fh Consensus here means full consensus, 
even without small extremist clusters. Further on, there is evidence that the 
consensus threshold for the HK model only remains to be (very roughly) about 
0.2 if the average number of links scales with the number of agents if these 
numbers go to infinity. Otherwise it switches to 0.5. In (37) the DW model was 
again set on a Barabasi- Albert graph and gives evidence that the number of final 
clusters scales now with the number of agents and not only with i. (See also 
([H) and |39l ) for results and some different versions of DW dynamics on scale- 
free networks.) In (j4p| ) a combination of multidimensional opinion dynamics on 
scale-free networks is done. 

The extremism version of the DW model proposed in (|24l ) has been set on 
small world graphs a la Duncan- Watts leading to the fact that the drift to one 
extreme only appears beyond a critical level of random rewiring of the regular 
networ k ( 1 4 if) . 

In (|23r ) a general model has been proposed which includes the DW and the 
HK model as special cases. A fixed number of m agents meet at each time step 
and perform a HK step. There is evidence that the transition from the HK to 
the DW model is somehow smooth. 

Finally, in (|36i ) it is emphasized that only a static network may modify the 
dynamics but also the choice of random communication partners in the agent- 
based DW model may be important. Such a choice is called communication 
regime. It is shown that enforcing as well as preventing consensus is possible 
for a for huge e-phase if it is possible to manipulate who meets whom at what 
time. 



Convergence parameter, cautiousness or quality The original agent- 
based DW model in (0) contained a convergence parameter \jl e]0, 0.5] which 
controls how far an agent moves to the opinion of the other if they are closer 
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than e to each other. The model presented here is the case fi = 0.5. For lower 
/i agents make shorter jumps towards the other opinion. Therefore, it can also 
be called cautiousness parameter. In a brief analysis of the basic model this 
parameter has been shown to affect only convergence time. Low [i leads to 
longer convergence time. But in it is shown by simulation that lower \x 
also reduces the sizes of the minor clusters. So more cautious agents are not 
that likely to produce minor clusters. Further on, /i can have an interplay with 
other parameters when they are introduced, e.g. the number of agents which 
meet in a time step in |23|) . It is possible to include the cautiousness parameter 
in a density based version of the model (36). In |4^) there is a study of the 
DW model on a growing scale-free network with quality differences of opinions. 
High quality is modeled as low n, which corresponds to high cautiousness. But 
quality is assigned to an opinion not to an agent. So, an agent's cautiousness 
changes when she changes her opinion. 

Miscellaneous In (@) the agent-based HK model is analyzed for biased 
confidence intervals. So, for instance agents look more to the left than to the 
right. This can lead to consensus at more extreme positions. 

The agent-based DW model has been extended two times to smooth bounds 
for the confidence interval, called the relative agreement model |24h and smooth 
bounded confidence model (|44l ) . The idea is that the attraction of agents should 
be lower if they are far from each other but still not totally zero. Results 
remained quite similar. Structural effects remain to be discovered. 

Another extension is the introduction of repulsive forces when agents are far 
away from each other. This is studied in the DW setting in (fiBT ) motivated by 
social judgment theory. There were absorbing boundaries at the extremes and 
in the center dynamics were similar to the usual clustering. 

The agent-based HK model has been extensively studied with other means 
than the arithmetic mean in (j46l ) leading to very similar clustering behavior 
except for what is called a random mean. But positions of clusters change a lot 
in comparison to the original model. Proofs of covergence are also possible for 
abstract means. 

Communication strategies of agents in the agent-based DW model were a 
topic in ((IB). Agents may be balancing and try to find communication partners 
from different directions or they may be curious and try to find communication 
partners of the same direction after a successfull communication. Both lead to 
better chances for consensus. But there is an interplay with the cautiousness 
parameter fi. For balancing agents it is better to be cautious to enhance chances 
for consensus, but for curious agents it is not good to be cautious. 



Urbig (|47l ) argues that it is important to distinguish opinions and attitudes 
(e.g. about a product). The attitude is a continuous variable which each agent 
holds internally. The attitude becomes an opinion when verbalized and com- 
municated. With this process it becomes discrete. In this terms this paper is 
about 'attitude dynamics'. 

From a more philosophical point of view the topic of finding the truth is 
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incorporated in the HK model in (|48l ) as an agent who automatically sticks to 
the 'true' value. Questions are then, how many agents who go for the truth do 
we need such that all agents end up finding the truth. In (|49l ) the idea is used 
for the DW model. 

5 Conclusions and open problems 

The models of continuous opinion dynamics under bounded confidence have 
been shown to capture some phenomena which have are not covered by tra- 
ditional binary opinion models. These phenomena are dynamic formation of 
different numbers of clusters with characteristic location and size patterns, and 
drifting of clusters when bounds of confidence are heterogeneous. Also from 
a modeling point of view continuous issues are different from discrete issues. 
Continuous opinions are more related to negotiation problems or fuzzy atti- 
tudes which do not really match with a yes or no decision. 

The master equation approach has lead to the density-based models which 
give better insight on attractive states of the agent-based models. The numer- 
ical integration of the master equation is something between a full analytical 
understanding and the analysis by Monte Carlo simulations. But on the other 
hand it is still important to match Monte Carlo runs with computations of den- 
sities because we are not only interested in dynamics in the thermodynamic 
limit but especially in dynamics of finite populations (of various sizes) . See also 
for results when finite population sizes are essential for opinion dynam- 
ics. Finite size effects play an important role even in very large populations 
(especially in the HK model) but on the other hand the density-based mod- 
els delivers attractive final cluster patterns for the agent-based model even for 
very low numbers of agents. But as ususal for lower population sizes a Monte 
Carlo run does not always match the attractive cluster pattern predicted by the 
density-based model. A quantification to what extend the density-based model 
matches with the agent-based models of a given population size would be a great 
improvement. 

From an analytical point of view the most important remaining issue is 
the understanding of the impact of the initial profile or initial density on the 
final outcome. The task is to classify initial profiles which lead to essentially 
different cluster patterns. The next important problem is to understand the 
drifting phenomena under heterogeneous bounds of confidence. Here, the master 
equation outlined above seems to be rewarding. But still the universal scale for 
bifurcation points is not proven today. Even the convergence to a stable cluster 
formation itself is not formally proven for the density-based models and the 
agent-based models under heterogeneous bounds of confidence. Trying to prove 
all that should also lead to more insight on dynamics. 

From a theoretical modeling point of view, it would be interesting to classify 
all models of continuous opinion dynamics which show stabilization of opinions 
in clusters. 

From an application point of view, a deeper analysis of continuous opinion 
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dynamics may lead to a better design of decision and discussion processes in 
large populations. Fostering the chances of finding consensus is a good starting 
point for research about extensions. The next level is then to foster consensus 
within the bounds of not too large convergence times. 
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